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We present a study of a gas of bosons confined in one dimension with Feshbach resonant interac- 
tions, at zero temperature. Unlike the gas of one dimensional bosons with non-resonant interactions, 
which is known to be equivalent to a system of interacting spinless fermions and can be described 
using the Luttinger liquid formalism, the resonant gas possesses novel features. Depending on its 
parameters, the gas can be in one of three possible regimes. In the simplest of those, it can still 
be described by the Luttinger liquid theory, but its Fermi momentum cannot be larger than a cer- 
tain cutoff momentum dependent on the details of the interactions. In the other two regimes, it 
is equivalent to a Luttinger liquid at low density only. At higher densities its excitation spectrum 
develops a minimum, similar to the roton minimum in helium, at momenta where the excitations 
are in resonance with the Fermi sea. As the density of the gas is increased further, the minimum 
dips below the Fermi energy, thus making the ground state unstable. At this point the standard 
ground state gets replaced by a more complicated one, where not only the states with momentum 
below the Fermi points, but also the ones with momentum close to that minimum, get filled, and 
the excitation spectrum develops several branches. We are unable so far to study this new regime 
in detail due to the lack of the appropriate formalism. 
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I. INTRODUCTION 

A gas of one dimensional interacting bosons is an in- 
teresting system to study, and we have in our possession 
a number of techniques to study it exactly. One of the 
first papers to study this system, Ref. demonstrated 
that with infinite repulsive interactions, it is equivalent 
to a gas of noninteracting spinless fermions. A subse- 
quent paper Ref. showed that even if the interactions 
are not infinite, as long as they can be approximated by 
short ranged delta function type interactions, the gas of 
interacting bosons can still be studied with the help of 
a method closely related to Bethe Ansatz. Its resulting 
behavior is in many ways equivalent to a gas of interact- 
ing fermions. In particular, two particles cannot have the 
same momentum when far away from each other, which is 
a sort of a Pauli principle for bosons. In the ground state 
all the momentum states below the Fermi momentum are 
occupied, while the rest are empty. 

Subsequently it was demonstrated in Ref. Q that the 
excitations of this gas around the Fermi points could be 
described by the effective "bosonized" theory which was 
later named Luttinger liquid theory. It is applicable to 
both systems of interacting bosons and fermions in one 
dimension. 

Interest in such systems has been revived recently 
due to the emergence of new experimental techniques 
where the one dimensional gas of interacting bosons can 
be manufactured in a laboratory with ultracold bosonic 
atoms as its constituent particles 4] . For the first time it 
has become possible to study these gases experimentally. 

A typical interaction between ultracold atoms is due to 
van der Waals forces, which are relatively short ranged 
compared with typical interatomic separation in a cold 
atom experiment, thus it is not impossible to conjecture 
that Bethe Ansatz techniques may be applicable to study 
the behavior of a gas of such atoms in one dimension. It 



has also become possible to control the interactions with 
the help of Feshbach resonances p| . Feshbach resonance 
allows to bring a pair of atoms in a resonance with a di- 
atomic molecule, whose binding energy can be controlled 
by the application of magnetic field. The binding en- 
ergy can be both negative, in which case the atoms will 
want to release their energy and form a stable molecule, 
or positive, in which case the molecular lifetime is finite 
and the molecule eventually decays into a pair of atoms. 

It is well known that, in one dimension, if the 
molecule's energy is negative, and if the atoms are 
bosons, then the system is unstable and collapses into 
a state whose energy is negative and is proportional to 
the square of the number of particles. Thus the energy 
per particle is essentially minus infinity in the thermo- 
dynamic limit. This situation cannot be experimentally 
realized. 

If, however, the energy of the molecule is positive, then 
this results in the new type of atomic gas not studied in 
the literature before. The study of the gas with resonant 
interactions is the subject of this paper. 

The Hamiltonian of such a gas can be written in the 
following form |j| 
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where a', b\ and a, b are bosonic atomic and molec- 
ular creation and annihilation operators. Here the in- 
teractions, given by the Feshbach term linear in b and 
quadratic in a, are automatically short ranged. The 
mass of atoms a is m a , while the molecules 6, being able 
to covert into two a-particles, must twice as large mass. 
Otherwise, Eq. (JTJ would violate Galilean invariance. eo 
is the parameter which controls the binding energy of the 
molecules. 
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We note that the Hamiltonian Eq. includes both 
the Feshbach g term and the delta-function interaction c 
term. Although often added when studying Feshbach res- 
onant gas in higher dimensions, the c-term can be shown 
to be unnecessary in three dimensional space [6|. It can 
be absorbed into the redefinition of eo and g. However, 
in one dimension the situation is different, and the delta- 
function term c is crucial. If the c term is omitted in 
Eq. QJ, then there exists a real bound state of atoms 
with negative binding energy at all values of eo and g. 
This particular situation does not arise in three dimen- 
sions. It is related to the fact that in one dimensional 
space (as well as in two dimensions) there exists a bound 
state in an arbitrary weak attractive potential. A bosonic 
gas with real bound state of atoms is unstable, as was 
discussed above. 

Fortunately real atoms always repel each other at short 
distances, and this is encoded by the c-term in Eq. IjTjI. 
This term is necessary to create a situation where real 
bound states are absent, while resonances, molecules with 
positive binding energy and finite lifetime in resonance 
with atomic gas, are present. 

A version of the Hamiltonian Eq. in higher than 
one dimensional space was analyzed in Ref. ;7j. It is well 
known though that the one dimensional physics is usually 
quite different from its multidimensional counterpart. In 
one dimensional space a similar Hamiltonian was studied 
in Ref. However, in that paper the a-particles were 
spin-1/2 fermions, and the analysis presented in that pa- 
per was based on the bosonization techniques, which do 
not capture the excitation spectrum beyond the vicinity 
of the Fermi points. In this paper we study Eq. ex- 
actly, and our techniques do not rely on bosonization. 
So we are able to obtain the exact excitation spectrum 
of Eq. 0. 

The question of whether the techniques presented here 
are also suitable to study the gas of spin-1/2 fermions 
with resonant interactions, as the one in Ref. |8|, is still 
open. The main difficulty seems to be the presence of the 
spin, which could complicate the Bethe Ansatz solution 
or even make it impossible to construct. We leave this 
question for future work. 

In this paper we construct solution to the problem de- 
fined in Eq. £Q| using Bethe Ansatz approach originated 
in Ref. 0. These authors study Eq. (Q with g = 0, so 
we extend their techniques to the g ^ case. It is im- 
portant to stay away from real bound states of atoms in 
our construction, to avoid the collapse of the gas. The 
condition for the absence of real bound states in Eq. 
(we will show it later in the paper) reads 
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> 0. 



(2) 



Thus we study Eq. supplemented by this condition. 

The bosonic Pauli principle for one dimensional inter- 
acting bosons of the sort studied in Ref. (that is, 
Eq. m with g — 0) demands that no two bosons can 
asymptotically occupy the same state. As a result, the 



bosonic system Eq. Q at g = can be described by the 
density of states p(A), where A is the asymptotic momen- 
tum (momentum when the particles are far away from 
each other). The ground state corresponds to all the 
states up to the Fermi momentum ±Xf filled, and the 
rest empty. The value of Xf depends on the total den- 
sity D of system (number of particles per unit length), 
with D(\p) being a monotonously growing function of 
its argument. The excitation spectrum e(A) can also be 
constructed 9] . It is a monotonously growing function of 
A at A > 0, negative at |A| < X F and positive at |A| > A^- 
It describes particle and hole excitations. In the vicinity 
of Xf it can be replaced by a linear function of A well 
described by the Luttinger liquid theory. 

We show that the Pauli principle also holds for the 
Feshbach Hamiltonian Eq. (JTJ when applied to the atoms 
|1C| . However, the excitation spectrum e(A) is more com- 
plicated than at g = 0. Three regimes of the resonant 
Hamiltonian can be identified, depending on the range of 
values of c. These are 



1. c< 





3. c> 




(3) 



In all three regimes, g is assumed to be in the range such 
that Eq. 101 holds. 

If the total density D is small, so that A^ is much 
smaller than the resonant momentum (momentum where 
an atom scattering off another atom with opposite and 
equal momentum would be in resonance with the di- 
atomic molecule), then the density of states p(X) and 
the excitation spectrum e(A) practically coincides with 
that at g = (delta function interacting Bose gas) stud- 
ied in Ref. Q and However, as the density of the 
gas is increased, the density of states and the excitation 
spectrum develop new features. 

In the simplest of the three regimes, where c ~ 

eo/m a , the excitation spectrum remains roughly sim- 
ilar to that of the delta function interacting Bose gas, 
even at higher density. However, as D is increased fur- 
ther, Xp can never exceed a certain cutoff momentum 
^cutoff p As Af approaches the cutoff momentum, the gas 
accommodates the higher density by increasing its den- 
sity of states p(A), as opposed to by increasing its Fermi 
momentum Xf- As a function of A, the density of states 
remains roughly constant in this regime. A™ can be 
calculated in a closed form, but the expression for it is 
not simple. It can be estimated to be A^ toff ~ m a c. 

In the regimes of small or large c, the behavior devi- 
ates even more from that of the delta function interact- 
ing gas. As the density is increased, e(A) ceases to be a 
monotonous function of A. Instead, it develops a local 
minimum above the Fermi point at a momentum such 
that the particles moving at this momentum are roughly 
at resonance with the particles in the Fermi sea. (see 
Fig. EH which depicts the case where c is small. The 
function plotted there is e(A) = e/X% where A = X/Xp-) 
This minimum cannot be captured by Luttinger liquid 
theory as it develops a finite distance above the Fermi 
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momentum and not in its immediate vicinity. 

As the density is increased further, the minimum in 
e(A) gets lower. At a certain critical density D CIlt the 
minimum touches zero, and above this density the mini- 
mum is now negative, Fig. 1121 At this point the ground 
state which resulted from filling the states with A between 
— Xf and Aj? becomes unstable and the states around the 
minimum also start getting filled. The actual value of the 
critical d ensity c an be found analytically only if c is small 
or c< y / eo/m^, and it is given by 



"crit — 



9m a ej] 
2 (ce - g 2 ) 



(4) 



The denominator of Eq. (J3J is positive due to Eq. J5J. 
The Fermi momentum of the critical gas is given by, again 
at small c 



X F = (6m 2 a (ce ~g 2 )y. 



(5) 



The case of large c is harder to study analytically, so 
we do not have analytic equivalents of Eqs. (@J and 10, 
but on the qualitative level the behavior of the excitation 
spectrum is similar in both cases. 

Unfortunately, we are unable to determine precise 
properties of the ground state and excitations once the 
system crosses over into this new high density regime. 
This is due to the fact that the equation which describes 
the new ground state is nonlinear, and it is not clear 
how to go about solving it, cither analytically or numer- 
ically. So in this paper we limit ourselves with specula- 
tions about the properties of this high density state. 

Which of the three regimes can be realized in the exper- 
iment depends on the properties of a particular Feshbach 
resonance. This is an important question, but we will 
not attempt to study it in this paper. 

Finally, we address briefly a possible criticism which 
can be brought against the technique we use in this pa- 
per. This technique, taking into account merely pair-wise 
interactions between the atoms, does not explicitly take 
into account that the atoms spend some part of their time 
bound into molecules, and that a molecule can interact 
with another atom in its vicinity. We observe, however, 
that the density of the molecules is very small, at least 
in the small and large c regimes and at total densities 
D not exceeding too much the critical density at which 
the phenomena discussed in this paper take place. This 
is due to the fact that few atoms are in resonance with 
each other in the regimes of interest here. Thus inter- 
actions between the molecules, in the regimes of interest 
here, can be neglected. We discuss this in more detail at 
the end of Section ITTll and in the Appendix [E] 

The rest of the paper is organized as follows. The 
Bethe Ansatz techniques require the knowledge of the 
two-body scattering matrix S. Section [H] is devoted to 
constructing the scattering theory in one dimension and 
applying it to Eq. on the two-body level. This sec- 
tion can be skipped at a first reading, and its only result 



important for further applications is the two-body scat- 
tering matrix Eq. I|27|) for the Hamiltonian Eq. 

In the next section lTTll we construct Bethe Ansatz many 
body wavefunctions out of the two-body scattering ma- 
trix S. The procedure is fairly standard and is included 
for completeness. 

The next section Hvl is devoted to the Bethe equations 
which allow to compute the density of states and the ex- 
citation spectrum of the Bose gas. Although standard, 
these equations involve the scattering phaseshift 5, re- 
lated to the S'-matrix as in S = e 2lS . The phase shift for 
resonantly interacting gas takes a new form which follow 
from Eq. 1|27|) . Thus in the remainder of that section we 
discuss qualitatively what behavior we can expect from 
5, identifying the physically distinct large c and small c 
limits, depicted on Figs. and The functions 9(A) 
shown on these figures is related to S by Eq. i|33|) . 

Armed with these preliminaries we proceed to con- 
struct the ground state and the excitations of Eq. 
in section It is in this section that the main results of 
this paper are derived and discussed. 

Finally, in section Ivll we discuss possible structure of 
the ground state and excitations at densities above the 
critical density. 

Finally, several appendices discuss technical details rel- 
evant for the material presented in this paper. 

One of the most detailed discussions in the literature 
of the standard one dimensional interacting Bose gas , 
without resonances, is presented in the book Ref. |ll| . 
Throughout this paper we extensively borrow their no- 
tations and use their results when needed. 



II. ONE DIMENSIONAL SCATTERING 

We begin this paper with constructing the two-body S- 
matrix for the Hamiltonian Eq. , which will be of use 
for the subsequent solution of the problem using Bethe 
Ansatz. Given a lack of discussions in the literature, it 
seems appropriate not only to construct it for this spe- 
cific Hamiltonian, but also to discuss its general features, 
which is the subject of this section of the paper. 

Quantum mechanics textbooks typically spend sub- 
stantial time discussing scattering in three dimensional 
space. At the same time, scattering in one dimension is 
discussed only briefly. One dimensional scattering is in 
many regards similar to the one in three dimensions, yet 
there are also crucial differences. 

Consider a particle in one dimension moving in a sym- 
metric potential U(x) — U(—x) which quickly vanishes 
at large distances, U(x) — s- if |x| tq, 



1 d 2 jj 
2m dx 2 



+ U(x)ip = Eip. 



(6) 



Here, as throughout this paper, the Planck constant H is 
set to 1. 

Typically for a problem defined in Eq. © we set up 
a transmission-reflection problem. It is defined as a so- 



4 



lution to Eq. © with the asymptotic behavior of ij>(x) 
given by 



ip(x) 



e lkx + re~ ikx , 



= te ikx , 



x < -r , 
x > r . 



(7) 



Here r and t are transmission and reflection amplitudes, 
with transmission and reflection probabilities given by 



T 



\t\ 2 , R 



However, for a symmetric potential, 



it is advantageous to set up a scattering problem defined 
in a way analogous to the scattering theory in 3D, 

4>{x) = e lkx + f s e- lkx - f A e- lkx , x « -ro, 
iP( x ) = e lkx + f S e lkx + f A e lkx , x » r . (8) 

Here fs and Ja are symmetric and antisymmetric scat- 
tering amplitudes (similar to the partial wave scatter- 
ing amplitudes in three dimensions). It is clear that 
r = fs — /a, and t = 1 + fs + f a- Another way to set 
up a scattering problem is via an S-matrix. It is defined 
via a solution to Eq. © with the asymptotic behavior 



tp(x) 



-ikx 



S s e 

+ s, 



— ikx 



se 



ikx 



x < -r , 
x > r 



for the scattering in the symmetric channel, and 

x < -r , 



Akx 



-ikx 



e" v " + Sag 

-ikx q „ikx 



S A e l 



x > r 



(9) 



(10) 



in the antisymmetric channel. By constructing ip(x) + 
ip(—x) and ip(x) — ip(—x) out of Eq. © it is easy to see 
that 



f — ^ s ~ 1 f — ^ A ~ 1 



(ii) 



By probability conservation the scattering matrix is uni- 
tary. As a result, the scattering matrix can be repre- 
sented in terms of phase shifts 



S s = e 2l5s , S A = e 2lSA - 



(12) 



Another convenient representation of the scattering am- 
plitude takes the form 



fs 



1 



ikF s (k 2 ) -1 



i 



(13) 



Here Fs{k 2 ) and Fa (A: 2 ) are real functions of its argu- 
ment, Taylor expandable in powers of k 2 at small k. The 
fact that Re(/g A ) — — 1 is simple to establish. It follows 
immediately from the unitarity of the S-matrix, or from 
|1 + 2fs,A\ 2 = 1- It is somewhat more difficult to estab- 
lish the form of Im(/^^), this is done in Appendix El 
All of these results closely parallel well known facts from 
the theory of 3 dimensional scattering (see, for example, 

G2)- 

In this paper, we are interested in scattering of bosons, 
which, due to the symmetry under particle exchange, do 
not scatter at all in the antisymmetric channel. So, we 



will not study scattering in the antisymmetric channel 
any further. To simplify notations, in what follows we 
will drop the subscripts S and A and use the notation 
such as / = fs, S = Ss, and F = Fs- 

If we are interested in scattering at sufficiently low en- 
ergy we can replace the function F(k 2 ) by its value at 
zero momentum, F(0) — —a. Here a is usually called 
the scattering length. Thus the low energy scattering 
amplitude and S-matrix take the form 



/ = -- 



1 



1 



S = 



iak — 1 



1 



(14) 



It is important to keep in mind that it is not enough to 
take I: « r 1 for this simplification to work. In fact, the 
condition specifying how small k must be for Eq. i|14|) 
to work depends on the details of the potential. In all 
the examples of the potentials of interest in this paper, 
k will be taken much smaller than ro, yet one will often 
have to go far below that scale for Eq. I|14f> to hold. Yet, 
for arbitrary potential there is always some scale so that 
below it Eq. ijTH is applicable [To|. 

Notice that in the limit k <C a -1 , the scattering am- 
plitude becomes / = — 1, and the S-matrix is S = — 1. 
This means that the transmission coefficient, t = 1 + /, 
vanishes at low energy, while the reflection coefficient 
r = / = — 1. Thus we arrived at a result of utmost 
importance in this paper: at very low energies the scat- 
tering in an arbitrary potential in one dimension results 
in total reflection (with one exception, see 01 )• Stated 
in this form, this is well known in q uantum mechanics lit- 
erature, see, for example, Ref. |12| . problem 5 in section 
25. 

The scattering phase shift is related to the function F 
via a simple formula 



cot<5 = -kF(k 2 ) w fee 



(15) 



where the last approximate equality holds at sufficiently 
low momenta. Thus at \ka\ <C 1, S = | if a is positive 
and 6 = — 77 if a is negative (we restrict S to the interval 
— § to J). Either way, this means that scattering at very 
low energies is always in the unitary limit. This result is 
specific to one dimension and has no equivalent in three 
dimensional scattering. 

An important potential to study is the delta-function 
potential U(x) = c 5(x), where c controls its strength. 
It is attractive at c < and repulsive at c > 0. Let 
us calculate scattering amplitude in this potential. It is 
easiest to do by first computing the phase shift. We look 
for a solution of the Schrodinger equation 

1 cPjijx) k 2 
~ 0^ a 2 + c 6(x)i>{x) = —ip(x) (16) 
2m ax^ 2 

in the form ip( x ) — cos(kx + 5) for x > and ^>{x) = 
cos(—kx + 5) for x < (this coincides with Eq. © up 
to an overall phase, and of course, in the case of delta- 
function, r =0). Direct substitution gives 



arctan ■ 



mc 



(17) 
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The scattering amplitude and S-matrix become 



/ = 



ik i 



s = 



ik + mc 
ik 



mc 



(18) 



Notice that at k = —imc the scattering matrix has a 
pole. This pole corresponds to the bound state with the 
energy E = k 2 /2m = —mc 2 /2 if c < and the potential 
is attractive. If, however, c > 0, then this pole does not 
correspond to the bound state and instead describes what 
is called in the literature a virtual bound state 0] ■ The 
existence of the bound state at c < 0, and its absence at 
c > 0, can also be verified by the direct substitution of 
k = —imc into Eq. 

The result expressed by Eq. ijTBl) is quite remarkable. 
It implies that the scattering amplitude in the delta- 
function potential coincides with the low energy asymp- 
totics of the scattering amplitude in arbitrary poten- 
tial, with the identification mc = —1/a. This is why 
when studying many body interacting systems with typ- 
ical energies low enough for the approximation leading 
to Eq. I|14f) to hold, it is sufficient to replace the realistic 
interactions by the delta-function interactions, which are 
often easier to analyze. And this is why exact results ob- 
tained in Ref . 0, E| for the delta- function interacting one 
dimensional Bose gas are applicable for the whole range 
of bosonic systems in one dimension regardless of their 
interactions. 

In this paper, however, we will be interested in po- 
tentials which support a long lived quasistationary state 
(resonance) at positive energy E = eo. A potential of 
this kind is depicted on Fig. ^ 

At energies far below eo, the scattering amplitude 
of this potential reduces to Eq. I|14f> . just like for any 
other potential. However, at energies of the order eo 
the scattering amplitude no longer takes the simple form 
Eq. (J2J. Indeed, the scattering amplitude must reflect 
in its pole structure the presence of the resonance at eo- 
That means, it must have a pole at the complex value of 
the energy 



2m 



e - iT/2, 



(19) 



where V is the inverse lifetime of the resonance 01 • 



U(x) 




Notice that the potential depicted on Fig. can be 
made arbitrarily narrow (and correspondingly deep, to 
keep eo at a finite value) , so eo can always be made much 
smaller than l/mr^. Thus momenta of the order of k ~ 
y/meo can be made much smaller than r^ 1 . Yet Eq. I|14|) 
may already not apply. Instead we have to go back to a 
more general result Eq. (|13|l . 

To work out a typical scattering amplitude for a po- 
tential such as the one depicted on Fig. ^ is not easy. 
Fortunately, we do not need to work with an explicit 
potential. Instead, we recall that we are interested in 
atoms interacting via a Feshbach resonance, described 
by Eq. . Two-body version of the Feshbach Hamilto- 
nian, sometimes referred to as a Fano-Anderson model 
[l3l | or a single level in the continuum is 



H 



1 

2m 



dx d x a) d x a- 



g [a\0)b + 6 + a(0)] + e & f 6 + c a f (0)a(0). (20) 



Notice that the numerical coefficients, present in Eq. (JJJ 
due to indistinguishability of the particles, disappeared. 
The Fano-Anderson Hamiltonian describes a particle, 
created by a), moving a long a one dimensional line 
which, upon hitting the point x — 0, may turn into 
another particle, created by b* . The amplitude of this 
process is given by g. The 6-particle does not move and 
has a fixed energy eo- This particle is unstable and af- 
ter some time will turn back into an a-particle which 
will then be free to move away from the origin. In ad- 
dition, we add a delta-function repulsion from the origin 
for the a-particle, with strength c. As discussed in the 
introduction, in one dimension the presence of the delta- 
function repulsion makes a crucial difference (eliminates 
"parasitic" bound states), as we will see below. 

It is not unreasonable to expect that this model de- 
scribes an effective potential with long lived resonance 
whose energy is approximately eo (assuming that eo > 0). 
The lifetime of the resonance is controlled by g. If eo < 0, 
then we expect that it describes a potential with a bound 
state. We will see below that this picture is more or less 
correct (although the energy of the resonance is eo at 
small c only. At large c, it is at ei, defined in Eq. J2HJ0- 

To understand the Fano-Anderson model deeper, let 
us compute the scattering amplitude of the a-particles. 
The scattering amplitude is given by the T-matrix, which 
is closely related to the scattering amplitude. In one 
dimension, this relation is given by 



f(k) = ~T{k). 



(21) 



T matrix itself, for the Hamiltonian Eq. (|20l) . can be 
computed diagrammatically. Leaving the details of the 
calculation to the Appendix we present the answer 



FIG. 1: A symmetric potential U(x) of the typical range ro 
with the quasistationary energy level at some energy e - 



/(*) 



ik 



2m 



eo 



' k 2 

\ 2/u 



(22) 
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First of all, we see that this scattering amplitude con- 
forms to a general form Eq. l|13[l a symmetric scattering 
amplitude in one dimension must assume. 

Second, we observe that in the absence of Feshbach 
term, g = 0, this scattering amplitude reduces to scatter- 
ing in the delta- function potential, Eq. (|18J) . as it should. 

Third, we would like to study this scattering ampli- 
tude in the absence of delta-function term, c = 0. The 
scattering amplitude reduces to a simpler form 

f c=0 (k)=(^(^-e )-l) \ (23) 
\mg z \2m J J 

This indeed describes scattering via a potential with 
a resonance close to eo (assuming it is positive) with 
the lifetime controlled by g. By analyzing the poles of 

Eq. (|23[1 it is easy to see that one of them, if j < e 2 / \frn, 
is given by 

*1 „ eo _ is a /Zl. (24) 
2m y V 2e 

From here and from Eq. (|19|1 we deduce that the life- 
time of the resonance is inversely proportional to g 2 . 

However, the scattering amplitude Eq. I|23(l also har- 
bors in itself a bound state at arbitrary values of eo- In- 
deed, if eo is sufficiently large, then the energy of the 
bound state is simply given by k 2 /2m « —mg 4 /2e 2 ] , 
where we balanced the term linear in k in Eq. (J13) with 
— 1, neglecting the cubic in k term. If eo is lowered 
until it becomes large negative, then this bound state 
smoothly crosses over into a more obvious bound state 
at k 2 /2m = eo, while the resonance which existed at pos- 
itive values of e simply disappears (mathematically, its 
pole moves off into an unphysical part of the complex 
energy plane). This can be seen by finding the corre- 
sponding poles exactly. 

The existence of the bound state of Eq. (121)1) at c = 
even when eo is positive and large can also be understood 
by noting that the a-particles, if their energy is below eo, 
can only spend a very short time in the 6-state. Their 
scattering will then be dominated by the second order 
perturbation theory where a particle jumps into a highly 
energetic 6-state with amplitude g, and then jumps back 
down to one of the a-states. Such process always leads 
to effective attraction between a-states and the origin. 
Unlike the three dimensional space, in one dimension an 
attraction, no matter how weak, leads to a bound state. 
That's the bound state we observe here. 

A system of one dimensional bosons with pair-wise 
short range interactions which support bound states 
is unstable and its energy per particle in the ther- 
modynamic limit actually becomes infinitely negative. 
(Ref. |2| ■ See also discussion in the next section.) That's 
why we would like to study potentials when there are no 
bound states. 

Fortunately, the full problem Eq. I|2U|I with additional 
repulsive interaction c > 0, whose scattering amplitude is 
given by Eq. (|22|l , does not have any bound states as long 



as the condition Eq. J2J) is fulfilled. Yet it does possess a 
resonance. Introducing notation 

ci = eo-— , 25 

c 

we rewrite the scattering amplitude as in 

m = (-p^-^) 1 - (26) 

\mc k z — Zmei J 

Here eo > e\ > 0. Indeed, at small c and g, this 
amplitude has a pole at k 2 /2m « eo — ig 2 y /, m/2eo, 
same as Eq. (I24|) . At large c and d, the pole is at 
k 2 /2m w ei — i(g 2 / c 2 )^2e\/m. In either of the two 
cases, the lifetime of the resonance is large. 

We take the scattering amplitude Eq. I|26[> as the basic 
resonant scattering amplitude of interest in this paper. 
For future reference, we also write down the S-matrix 
corresponding to the amplitude Eq. I|2()l) 

ik k 2 — 2mep ■ ^ 

S(i)=2/(t) + l= ;; : g (27) 

me k 2 — 2mei 

The two limits of small or large c and g, where the 
problem possesses a narrow resonance, have an interest- 
ing interpretation. In one of the regimes, where c <C 
yeo/m, we observe that f(k) is very small almost every- 
where except in the vicinity of k — and k 2 /2m = e , 
where / w — 1. Recalling that the transmission ampli- 
tude t = 1 + /, we find that this regime corresponds to 
perfect transmission at all values of momentum except at 
very low momenta and at the resonance, where we have 
total reflection. This is the regime 1 of Eq. 10 • 

The other regime, c 3> y/eo/m, the regime 3 of 
Eq. J2J), leads to f(k) = —I almost everywhere except 
at k 2 /2m — ei, where / = 0. Thus it describes a system 
with total reflection at all values of energy except at reso- 
nance (which is at ei), where we have total transmission. 
It is not unreasonable to expect that this situation de- 
scribes the transmission/reflection amplitude in case of 
the potential depicted on Fig. ^ (with ei playing the role 
of eo on the figure). 

Whether any of the two limits is appropriate to de- 
scribe experimentally observed Feshbach resonances de- 
pends on the properties of a particular resonance. This 
question deserves further study, but will not be discussed 
in this paper. 

III. BETHE ANSATZ SOLUTION TO THE 
MANY-BODY SCHRODINGER EQUATION 

Consider a gas of bosons moving on a line and inter- 
acting via a pair-wise interaction U{x) 

1 N f) 2 I 

~2^d^ + £ U( Xj -x k )=E^. (28) 

j 3 N>j>k>l 
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The wave function if)(xx, X2, • ■ ■) must be symmetric un- 
der interchange of any of the coordinates. For simplicity, 
we set the mass of bosons to be equal to m a — 1; hence 
the coefficient 1/2 in front of the second derivative. Note 
that the reduced mass of a pair of bosons is now set at 

to = 1/2, 

so all the scattering formulae of the previous chapter 
should be understood as having this mass. 

Generally speaking, for a large number of particles N 
and arbitrary potential U(x) it is not possible to solve 
this equation. However, if U(x) is very short ranged, 
then the probability that three particles simultaneously 
interact with one another is very low, and we can as- 
sume that only two particles interact with each other at 
a given time. In this case, we can use Bethe Ansatz to 
construct many body wave function. In what follows, we 
will closely follow Ref. m a Pplymg Bethe Ansatz to 
this problem, emphasizing the differences whenever they 



ip(x 1 ,X2, ■ ■ ■) 

Tp(xi,X 2l ■ ■ •) 
4>(xx,X2, . . .) 



Here 5(A) is the scattering matrix in the potential U(x), 
computed by solving the Schrodinger equation Eq. © 
with to being the reduced mass of the particles, m = 1/2. 
To write Eq. I|29|l down, we add a factor of S((Xj — Afe)/2) 
for every elementary permutation between Xj and Xk (j < 
k) one needs to do to rearrange the list X\, X2, £3, • • • in 
the order of increasing x, Xj 1 < Xj 2 < . . .. 

Then we symmetrize the function ip to arrive at the 
symmetric wave function 

i> = '^2i>(xp 1 ,xp 2 ,...), (30) 
p 

where Pj represents permutations of 1, 2 . . . ,N. ip satis- 
fies Eq. (|28|) with the energy 

j 

This is proven in Appendix O This wave function is 
a generalization of a standard Bethe Ansatz construc- 
tion for the many body problem Eq. (|28|) with the delta- 
function potential [TT| . 

One notable property of this wave function lies in the 
following: if any two A happen to coincide, tfi = 0. In- 
deed, suppose Xj — Xk ■ Then it is straightforward to see 



are present. 

Usually one constructs these wave functions for the 
delta function interactions |ll|. However, as we saw in 
the previous section, there are other potentials which 
are, on the one hand, short ranged, and on the other 
hand, have interesting structure. The potential depicted 
on Fig. ^ is one example if it is made sufficiently short 
ranged (tq <C 1/k and also much less than the typical 
spacing between the bosons). Even better example is 
given by the particles interacting via Fcshbach resonance 
Eq. QJ. J us t as m Eq. I|28l) . in what follows we take the 
mass of the bosons m a in Eq. Q to be 1, which means 
that the mass to in the two-body scattering matrix 5 is 
now equal to the reduced mass of the bosons, or 1/2. 

Bethe Ansatz for the many body wave function which 
solves Eq. (I28f) or Eq. can be written down in the 
following manner. First we construct an unsymmetrized 
wave function i\) 



(29) 



that ip(x\,X2, ■ ■ ■) is antisymmetric under interchanging 
Xj and Xk- Two crucial properties of the S-matrix nec- 
essary to prove that are 5(0) = —1 and 5(A)5(— A) = 1, 
and the proof is discussed in more detail in Appendix IU1 
Thus ip, which is a symmetrized version of ip, vanishes. 
Therefore, all Xk must be different for the solution to the 
Schrodinger equation to exist. 

This property of the many body one-dimensional 
Schrodinger equations for bosons is sometimes referred 
to as a bosonic "Pauli principle" . Interacting bosons in 
one dimension must all have different momenta, as if they 
are fermions. While discussed in textbooks for delta- 
function potential |ll|. this property holds for arbitrary 
short ranged potentials (with one exception [T(ij|'). 

Finally let us discuss the importance of the 3-body 
processes neglected here. We observe that in the process 
of checking that the wave function Eq. (|30|) satisfies the 
Shrodinger equation, we neglected the following 3-body 
processes. Suppose two of the coordinates, say Xj and 
Xk, are taken to be equal to each other. And suppose a 
third coordinate xi is passed through this point. Then 
Eq. H30|) may violate the Schrodinger equation in this 
range of coordinates. 

When the interactions are given by delta function, then 
the integrability protects the 3-body processes from vi- 



3 »E 3 - x j x j^ xx < x 2 < ■ ■ ■ < x N , 



Ai — X 2 



X2 < X\ < X3 < . . . < X]y, 



e iT. } ^ S { S ( I , X2 < x 3 < Xl < x 4 < 



< x N , 
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dating the Schrodinger equation. More generally, if the 
interactions are short ranged, then the weight of the wave 
function when two points are within the range of the in- 
teractions To is negligibly small, and 3-body processes can 
be neglected. However, in our problem, due to the pres- 
ence of the resonance, we can imagine that the wave func- 
tion has a spike every time two particles in resonance ap- 
proach each other. Thus a substantial part of the weight 
of the wave function can be at points where Xj = Xk- 

Although one can imagine this may invalidate our so- 
lution Eq. (|3UII when the density D of the gas is so high 
that there is a finite density of the molecules in the sys- 
tem, we show in the Appendix [E] that in the interesting 
regimes of small c or large c, at the density of the gas 
around the point of instability (given by Eq. at small 
c, for example), the density of the molecules is still suf- 
ficiently small so that we can neglect 3-body processes 
and rely on Eq. (|30|l for constructing solutions to our 
problem. 

As the density of the gas is increased, the molecular 
density increases as well, and at some point the tech- 
niques discussed in this paper breaks down. However, 
this happens at densities well above the ones where the 
phenomena discussed here take place. 

Ideally it would be interesting to find a model, similar 
in its physical properties to Eq. which would be inte- 
grable so that the integrability would protect the Bethe 
Ansatz at all densities. Whether this is possible is not 
known. 



IV. BETHE EQUATIONS 

Not all sets of Aj are appropriate for the wave function 
Eq. (|3*U|l . The wave function must satisfy appropriate 
boundary conditions. As always, the specific form of the 
boundary conditions is irrelevant while it is easiest to 
study periodic boundary conditions. Note that if ip sat- 
isfies periodic boundary conditions, then tj} also satisfies 
them. Imposing periodic boundary conditions on -0, we 
find 



Aj — Afe 



= 27m,,-, 



(32) 



where hj are arbitrary integer numbers, and L is the 
size of the system. It is customary in the Bethe Ansatz 
literature to define a modified scattering phase shift 0, 
related to the standard phase shift 2iS = logS 1 in the 
following manner 



0(A) = 25(A) +7T, A > 0, 
6»(A) = -0(-A), A < 0. 



(33) 



The scattering length of the S-matrix of interest to us 
Eq. Q26[l. under condition Eq. J5J), is negative (this is 
related to the fact that the potentials of interest to us 
have no bound states). For negative scattering length, 
lim,x— vo (5(A) = — 7r/2 (see discussion after Eq. <|15[l . Thus, 



we find that 9(0) = and 9 is an antisymmetric contin- 
uous function. In terms of 9, Eq. Il.'>2[l becomes 



N 

/c=i 



Aj - Aj 



= 27m 



(34) 



where rij are integers if N is odd, and rij are half integers 
if TV is even. These equations are called Bethe equations 
in the literature. They provide restriction on possible 
values Xj may take. We need to solve them to find the 
energy spectrum of our problem. 

Let us recall several useful properties of Bethe equa- 
tions. It is shown in Ref. [llj that all the solutions Aj of 
the equations Eq. (|34|l are real provided that 

\S(X)\ < 1 when Im A > 0, 
\S(X)\ > 1 when Im A < (35) 

Indeed, Eq. l|T7Jl . together with Eq. J2J), satisfies these 
conditions. This can be shown with some straightforward 
algebra. Thus, in the absence of bound states (condition 
expressed by Eq. (J2J)), all solution to the Bethe equations 
are real. 

Second important property stems from the fact that 



dX 



> 0. 



(36) 



This is also true curtesy of Eq. @ , as can be checked di- 
rectly by examining the explicit form of the derivative of 
9, given by Eq. H4()|) . It is shown in Ref. that as a con- 
sequence of positivity of derivative of 9, the solutions to 
the Bethe equations Aj can be uniquely parametrized by 
the set of integers (half- integers) rij. Conversely, in the 
presence of bound states (violation of condition Eq. J2J) ) , 
the function 9 is no longer monotonous, so there could 
be multiple solutions to Bethe equations. 

Finally, it is shown in Ref. pl| that under the same 
condition Eq. I|3tj|) if rij > n^, then Aj > A^. If rij 
then Aj = A&. Recall that if Aj 



At, then the wave 



function tp vanishes. It follows that no two rij may be 
the same. 

To summarize, a set of distinct integers (half- 
integers) rij , under conditions of interest to us, uniquely 
parametrizes the allowed momenta Aj of the particles, 
and larger rij correspond to larger Aj. 

In the remainder of this section let us analyze the be- 
havior of 9 of use in subsequent discussions in this paper. 

For a delta function potential, the phase shift is given 
by Eq. H17fl . The phase 9, constructed with the help 
of Eq. is depicted on Fig. |21 It is a monotonously 
increasing function, reaching tt at +oo and — tt at — oo. 

For a resonant scattering amplitude Eq. I|26|l. the form 
of the phase 9(X) depends on the two possible regimes, 
discussed in text after Eq. (|27|l . In the regime of small c, 
it is schematically depicted on Fig. EI At |A| < Ao, where 
a convenient notation is introduced to be used extensively 
in what follows 



Aq = V2meo, Ai = \ / 2mei = \ j 2m ( eo J , (37) 
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FIG. 2: A plot of 9(X) for a delta function potential. 
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FIG. 3: A plot of 9(X) for a resonant potential with scattering 
amplitude given by Eq. I26L and with c small. Here An = 
V2meo- 



it closely follows the graph for 0(A) for the delta func- 
tion potential. However, at larger |A| the phase jumps 
by another 2tt, characteristic of resonances in scatter- 
ing theory (the phase shift 5 jumps by ir). Notice that 
0(±Ao) = ±27r. This describes a situation where at al- 
most all the energies 9 — it or 9 — 3tt, which corresponds 
to the total transmission. In the vicinity of A = and 
A = ±A , the phase goes though 9 — or 9 — ±2tt, 
which is total reflection. The total change in the phase 
0(co) - 6>(-oo) = 6?r. 

If c is large, then the phase 9 is depicted on Fig.0] The 
phase is either or 2tt almost everywhere except in the 
vicinity of A = ±Ai, Ai = v / 2mej", where 9 goes through 
7T. This corresponds to total reflection almost everywhere 
except at energy equal to t\. We note that at very large 
A, not depicted on Fig.0] 8 subsequently increases to 3tt. 
However, it does it very slowly and over a long interval 
of momenta (in the limit of very large c). We will not 
consider behavior of 9 at such large momenta. So for our 
purposes 9(oo) — 9(—oo) = An in this case. 

Finally, we repeat that to conform with Eqs. (|28|l and 
(Q, to in the preceding paragraphs should be chosen as 
the reduced mass of the particles there, or m = 1/2. 

In what follows, we will need the function 



(38) 



For delta function potential, this function can be easily 
computed, by differentiating Eq. (|17|l . The result is given 



FIG. 4: A plot of 6(X) for a resonant potential with scattering 
amplitude given by Eq. 1261 and with c large. Here Ai = 
\/2mei. 



by 



K(X) 



2c 



A 2 ' 



(39) 



For the resonant amplitude Eq. (|26() this function can 
also be computed directly, with the result 



K(X) = 2c 



A 4 + 4 (Aq - 3A 2 ) A 2 + 16A§A 2 
(A 3 - 4AA 2 ) 2 + c 2 (A 2 



4A?) 



(40) 



This formula is not very illuminating, except if used for 
numerical calculations. Instead, we will concentrate on 
the cases of small or large c. We notice that K, being 
the derivative of 9 depicted on Fig. [3] and Fig.0] is zero 
almost everywhere except in the vicinity of 2Ao and 2Ai 
respectively. Thus, it can be well approximated by 



K{\) 



2ci 



2c 2 



2c 2 



A 2 



c 2 + (A - 2A r 



(A 



2A r 

(41) 



for small c, with 



Cl 



A 2 

An ' 



' 2A 2 



9, which corresponds to Eq. I|41|l . is depicted on Fig. 
For large c, the approximation reads 

2c 3 2c 3 



K{\) = 



c 2 + (A — 2MY 



+ (A + 2Ar) 2 



(42) 



where 



c 3 = 2 



\2 \2 

A Q - A 1 



and the corresponding 9 is depicted on Fig. 0] 

If c is neither small nor large, a surprisingly good ap- 
proximation amounts to replacing K{\) by its value of 
A = 0, or 



9{\) 



4A§. 



(43) 



This approximation, although crude and not extendable 
to large A (since 9(\) < 3ir), nevertheless allows to cap- 
ture most features of the moderate c regime rather well, 
as we will see below. 
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V. GROUND STATE AND EXCITATIONS OF 
THE RESONANT BOSE GAS 

To construct the ground state of the system of N par- 
ticles governed by the pair-wise S-matrix S, we need 
to choose a set of Xj, which satisfies Bethe equations 
Eq. 1(34(1 while minimizing the energy Eq. (|H 1|) . with all 
Xj being different. It seems reasonable that for that we 
need to choose rij as small as possible. Together with the 
requirement that no nj are the same, the ansatz for the 
ground state reads 



N-l N-l 



N-l 



(44) 



This ansatz indeed works for the delta function poten- 
tial ^3 • & ^ e shown below that in case of resonant 
scattering at small or large c, this ansatz breaks down as 
soon as the density of the gas exceeds a certain critical 
value. Nevertheless, we would like to work out the con- 
sequences of Eq. <(44ll here, in the process obtaining the 
criteria when Eq. (I44J) holds. 

Following Ref. |llj |. we introduce function X(x) such 
that X(rij) — Xj, which satisfies, buy virtue of Eq. (|51Jl . 
the equation 



jv 



LX{x) + Y J ° 



fe=l 



X(x) - A fe 



2nx. 



(45) 



Differentiating this equation with respect to A, and in- 
troducing a function 



P(A) 



1 dx 
LdA' 



(46) 



where L is the length of the system, we find in the ther- 
modynamic limit N — > co 



Zir 



dpK{X-p)p{p) = ^, (47) 



where K was defined in Eq_. Eq. ljT7l) is called the 

Lieb-Liniger equation 0, [lj, and K can be referred to 
as the kernel of this integral equation. Here Xp = Xn = 
— Ai , is the highest momentum the particles in the ground 
state have. In this way, it is akin to Fermi momentum, 
thus the notation. p(p) has the meaning of density of 
Xj, with p(Xj) — l/L(Aj+i — Xj) in the thermodynamic 
limit. In other words, it is the density of states per unit 
momentum. The total gas density can be expressed in 
terms of p as follows 



D=» = 
L 



dp p(p). 



(48) 



By construction, p(n) must be a positive function. It 
is shown in Ref. [ll| that this implies that 

5(X - p) - ^K(X - p) 



must be a positive definite operator, with its smallest 
eigenvalue no smaller than l/(27rp max ), where p max is 
the maximum value of p(A) on the interval between —Xp 
and Af- 

A necessary condition for this operator to be positive 
definite is 



dX 



dp 



S(X-p)-^K(X-p) >0 (49) 



This is equivalent to 
1 

a7 



dX ffi < 1. (50) 



For the delta function potential whose phase 9 is depicted 
on Fig. [2 Eq. 1(50(1 is always fulfilled no matter what Xp. 
Thus, Xp can be arbitrary. And indeed, it is possible to 
show that Xf depends on D, and the bigger the total 
density D is, the bigger Xp is. 

However, the situation is drastically different for the 
resonant 6, depicted on Figs. and 21 The scattering 
phase 9 from Fig. [3] results in the relation 



Xp < An 



(51) 



At the same time, the scattering phase from Fig.0]results 
in 

X F < 2Ai. (52) 
In case of moderate c, we can use Eq. (|4*5|) to find 

A^<^|, (53) 

from which the estimate quoted in the Introduction fol- 
lows. As a result, the Fermi momentum of the resonantly 
interacting Bose gas in one dimension cannot be arbitrar- 
ily large. Increasing the density of the gas D increases 
the Fermi momentum A^. However, as A^ approaches its 
limit Eq. (JSTJ, Eq. or Eq. (JHSJ) the density of states 
p{X) starts increasing very fast, and D diverges to infinity 
as A approaches its upper limit. This is confirmed by nu- 
merical solution of Eq. 1)47(1. We will see below, however, 
that only in case of moderate c this limit can be reached. 
In case of small or large c, an instability develops in the 
ground state before the limit is reached. 

The excitations above the ground state governed by 
the Lieb-Liniger equation Eq. I|47() are controlled by the 
function e(A), which obeys 

e(A) - y f F dp K(X - p)e(p) = y - h. (54) 

Here h is the chemical potential and A 2 /2— h is the energy 
spectrum in the absence of interactions. This equation 
must be solved with the condition e(±Xp) = 0. This 
condition fixes h if density D is given. 
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Once the function e(A) is known, the energy of an ex- 
citation is given by e, while the momentum of this exci- 
tation is given by k, where 



k = A 



dfi 



A — fi 



(55) 



Positive e correspond to a particle, while negative e to a 
hole. 

It is proven in the literature (see Ref. that in 

case where K(X) is a monotonously decreasing function 
of A for positive A, then the solution to Eq. I|54|l is a 
monotonously increasing function of A for positive A. An 
example of such K(X) is given by Eq. I|39l) . or by the 
delta function potential. Then e(A) < when |A| < Xf 
and e(A) > when |A| > Ajr. It is consistent with the 
ground state where all the states with |A| < Xf are filled, 
and the rest are empty. The excitations are then particle- 
like for |A| > Xp or hole-like if |A| < Xf- 

The situation changes drastically if K(X) is no longer 
monotonous, as in the resonant case Eq. (|41|) or Eq. (|42[) . 
Then e(A) is no longer necessarily a monotonous function 
itself. 

Eq. I|47l) and Eq. I|54[l are integral equations which in 
general are hard to solve analytically. In case of delta 
function interactions, where K is given by Eq. (JSHJ), it 
can only be solved in limits of very large c or very small 
c. The corresponding solutions are described in Ap- 
pendix [D] 

In case of resonant interactions, where K is given by 
Eq. (|40|l . we are able to find the analytic solution in two 
cases: the case of moderate c where we can approximate 
the phaseshift by Eq. I|43|) . and in the case of small c only, 
where K reduces to Eq. I|41|l . We proceed to describe the 
solution in these cases. 



Ground State and Excitation of Resonant Gas: 
the Case of Moderate c 



When c ~ ^/to/n^a, it is possible to obtain the quali- 
tative estimate of the solution by employing the approx- 
imation Eq. H43fl . or 



K{X) 



2A| 
cX\' 



(56) 



Since K is a constant, solving the corresponding integral 
equations is straightforward. Substitution into the Lieb- 
Linigcr equation Eq. I|47|) gives 



1 - 



2A^A 



A F 



X\kc 



(57) 



Thus the density of states is a constant, and Xf cannot 
exceed its limit given by Eq. (|53l . We will see below, 
by solving the Lieb-Liniger equation numerically, that 
although a crude approximation, Eq. I|57|l captures the 
main features of the density of states: it is constant and 



it diverges to infinity as Xf approaches its upper limit. 
The total density is approximately 



7T 



1 



1X\Xp 



X\kc 



(58) 



It is equally easy to find the excitation spectrum by 
solving Eq. . We find 



e(A) 



2 2 ' 

where the chemical potential h is equal to 



Ap I '_ 4AqAj? 



cttX\ 



(59) 



(60) 



Thus the spectrum remains quadratic, as for the nonin- 
teracting Fermi gas. We can now find the Fermi velocity 
vf by using 



de 
dk 



x=\ f 



D 

~dTT 



X, 



1 



2AqA^ 
X\ttc 



(61) 



(see Ref. for derivation). Here k is the physical mo- 
mentum of the excitations, given by Eq. (|55|) . We can 
also find the Luttinger parameter 



.9 



dD 
r dXF~ 



1 



2Aq Xf 



Xfirc 



(62) 



The Fermi velocity and the Luttinger parameter both 
diverge as Xf approaches its upper limit. 



B. Ground State and Excitation of Resonant Gas: 
the Case of Small c 

When c is small, we can do far better than when c is 
moderate, and can obtain solution which is not just an 
estimate, but is correct in the limit of vanishingly small 

c. 

As we will see in a second, in case when c is small, we 
only need to consider such Ao that Ao A^ . Under this 
condition, for the interval |A| < A^ , we can neglect the 
second and the third term of Eq. I|41|) and the function 
K reduces to the one for the delta function interacting 
potential, given by Eq. I)39fl . with the substitution c\ — > 
c. The corresponding Lieb-Liniger equation, given by 
Eq. 1D1I) . has been studied in the literature for a long 
time. In particular, its solution for c being small can be 
found analytically. The solution is presented in Appendix 
IdI and the answer for p(X) in the interval —Xp < X < Xp 
is given by Eq. (|D3|) . 

Now if we would like to know p(X) beyond this interval, 
we can no longer neglect the second and third term of 
Eq. I|41l) . However, we can use the Lieb-Liniger equation 
to find p(A) for any A by using 



P(A) = 



Xi 



An 2 a 



d M *T(A-/i)Jl-| r + ^. ,(,:!) 



1 

2n' 
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Taking into account that c\ and c 2 are small, we can 
approximate if by a sum of three delta functions, 



K = (5(A) + 5{X - 2A ) + S(X + 2A ) 



to find 



p(X) = 



and 



2?rci 



1 - ^2-, -Xf < A < A F , 



A, 



'1 - 



"< A > " 2« 2 » 
2Aq — Xp < A <• 2Aq 



(A-2A )^ 



x 2 

F 



(64) 



(65) 



(66) 



numerically, that they give a rather good approximation 
to the actual solution. 

Looking at the expressions for e(A) we indeed find that 
it is not a monotonous function of A. In fact, it has a min- 
imum approximately at A w ±2Ao. The minimum occurs 
around momenta whose particles are in resonance with 
Fermi sea. This is because two particles with momenta 
2Ao and have a relative momentum Ao in their center 
of mass frame. Thus they are in resonance. 

Moreover, it is easy to convince oneself that if Xf is 
increased at a fixed A , at a certain critical value of Xp 
the minimum becomes negative and the excitation spec- 
trum becomes unstable. The critical value of the Fermi 
momentum when the minimum crosses zero can be also 
approximately found by equating Eq. H70|) with zero and 
is given by 



A, 



(A + 2A r 



'1- 



P(A) = _ 

2-7TC2 

— 2Aq — Xf < A < — 2Aq + Xp- 



(67) 



Outside these three intervals within small c approxima- 
tion p(A) = l/(27r). The total density of the condensate 
coincides with its delta function value, or with Eq. I|D4|) 
with c being replaced by ci , or 



4 Cl 



(68) 



Now we use the same technique to find e(A). Within 
the interval — Xp < X < Xp the solution is simply given 
by the corresponding delta-function interaction solution, 
Eq. (|D7(1 . Outside this interval, we can use the approxi- 
mation Eq. (|64|) . Thus we find 



, |A|<Ai 



(69) 



A, 



[i2 Cl x 2 y 



(73) 



Upon substituting for c\ and expressing Ai, Ao in terms 
of the parameters of the Hamiltonian, we find Eq. JSJ). 
Notice that this critical value of Xf is much less than 
Ao, thus the approximation used to derive these results 
indeed holds. This leads to the critical value of density 



9A|^ 3 
4ci 



(74) 



This is equivalent to Eq. |0J. Above this critical density 
the Bose gas with resonant interactions as in Eq. Q with 
c small no longer has a simple ground state described in 
this section of the paper, as the function e(A) becomes 
negative at momenta higher than Xf- 

To check these results, we are going to compare them 
with the numerical solution of the Lieb-Liniger and exci- 
tation spectrum equations. 



Numerical Solution 



and 




X 2 X 2 X 3 
e (A) = L - ±£L _ ±£. 
' 2 4 6ci 



(A - 2Ao) z 

X 2 
F 



1 - 



(A + 2A r 

X 2 
F 



— 2Ao — A F < A < — 2Ao + Xp. 

Outside either of the three intervals 

A 2 



e(A) = 



(70) 



(71) 



(72) 



We emphasize that these formulae break down around 
the ends of the three appropriate intervals. We will see in 
the next section where we solve the Lieb-Liniger equation 



For the purpose of numerics, it is advantageous to pass 
to dimensionless variables by the appropriate rescaling 

A = Ai?A, /i = Xpfl, c = Xpc, e = Xpt, h — Xph. (75) 

p does not need to be rescaled. Then the Lieb-Liniger 
and the energy spectrum equations become 



1 



2tt 



1 



P(A) - — / dp K(X - p)p(p) = — 



2tt' 



and 



<(A)- j- I ,lf, A'(A-/W) = ^--/< 



(76) 



(77) 



Here K denotes the functions given in Eq. (|41|l or 
Eq. (|4*2*|l . but with c\ = ci/Xf, c 2 — c^IXf, £3 = c^/Xf, 
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FIG. 5: p(A) as a function of A. Here K is taken from Eq. 11 1 ( )ti 
with Ao = 2, Ai = 1, and 5=5. Dashed curve is the analytic 
estimate Eq. ^57|. 
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FIG. 6: e(A) as a function of A. Here K is taken from Eq. 14011 
with Ao = 2, Ai = 1, and 5=5. Dashed curve is the analytic 
estimate Eq. 159H . 



Ao = Ao/Af , and Ai = Xi/Xf replacing the variables 
ci, c 2 , c 3 , Ao, Ai in them. 

These dimensionless equations have been solved nu- 
merically. We used a simple discretization procedure 
where the interval between —1 and 1 was divided into 
1200 intervals, and the integrals in the equations were 
computed via a straightforward Newton discretization. 
The resulting system of linear equations was solved by 
inverting the kernel matrix. 



1. Numerical solution for moderate c 

We take Ao = 2, Ai = 1, c = 5. Under these conditions, 
the upper limit Eq. I|53|) for Xp is around 2. We solve the 
integral equations Eq. (|47|l and Eq. I|54|l for Xp = 1.96, 
close to its upper limit. Fig.[5]shows the density of states 
(solid line) and the analytic estimate Eq. I|57l) (dashed 
line). Although not coinciding, the estimate captures 
the fact that p is close to a constant, and that it diverges 
as Ai? approaches its upper limit. 

Fig. H3 shows the excitation spectrum, both numerical 
and the analytic estimate Eq. (|59|l (shown as a dashed 
line) . 
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FIG. 7: p(A) as a function of A. Here K is taken from Eq. I|41[l 
with ci = C2 = 0.01, and Ao = 10. Solid line represents 
numerics, while the dashed line is the analytic prediction. 




FIG. 8: e(A). Here K is taken from Eq. (UTI with ci = 5 2 = 
0.01, and Ao = 10. The solid curve represents numerics, and 
the dashed line is the analytic prediction. 

2. Numerical solution for small c 

We take K from Eq. i|41[l . First, we consider the case 
when Ao » 1, or in other words the energy of the reso- 
nance Ao is much bigger than the Fermi momentum Xf- 

Fig. \7\ presents the graph of p(X) obtained numerically 
for the case when Ao = 10 and c\ = C2 = 0.01. For com- 
parison, the analytic prediction Eq. Q65[) is plotted as a 
dashed line which generally lies slightly below the nu- 
merical line. This graph is indeed very close to Eq. 16511 , 
except at A = ±1 where, as discussed in the Appendix ID1 
the analytic solution breaks down. 

Fig. OH presents the graph of e(A). The analytic pre- 
dictions, which follow from Eq. (|M9"|l and Eq. l(72l . also 
work rather well. In fact, for |A| > A^ , we used Eq. i|72|) 
as opposed to a more precise Eq. 1D8|) . Had we used 
this latter expression, the analytic curve would have been 
completely indistinguishable from the numerical line. 

Finally, Fig. [5] presents the graph of e(k), with k being 
the rescaled momentum of the excitations, k — k/Xp, 
where k is given by Eq. (|55(l . For simplicity, we approxi- 
mated 5 in Eq. fB^f by S ~ Trsign (A) + Trsign (A — 2Ao) + 
7rsign (A + 2Ao). This graph also closely follows the ana- 
lytic prediction given by Eqs. CO, EHJ, and (El. We 
would like to remark that this graph shows that the ex- 
citation spectrum at A ~ ±Xp undergoes a sharp change 
in the slope. It is interesting that the Luttinger liquid 
theory concentrates on a tiny interval around the Fermi 
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FIG. 9: e(fc) as a function of the momentum of the excitations 
k, computed numerically. Here K is taken from Eq. lilt with 
ci = £2 = 0.01, and Ao = 10. 
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FIG. 10: e(A) as a function of A. Here K is taken from 
Eq. (gljl with ci = c 2 = 0.01, and A = 3.3. 



points where the spectrum can be linearized. It is clear 
that this theory cannot tell us anything about the spec- 
trum away from these points (the interval where Lut- 
tinger liquid theory applies is of course so small because 
of the small value of c). 

We expect that as the density of the gas is increased, 
or in other words, as Ao is decreased, the excitation spec- 
trum will develop minima. At the critical value of Ao, 
given by Eq. I|7'ljl. the minimum of the spectrum will dip 
below zero. At c\ = 0.01, this value is Ao_~ 2.89. We 
solve the integral equations numerically at Ao = 3.3. 

The density of states in the interval — 1 < A < 1 re- 
mains basically unchanged, as predicted. Let us now look 
at the excitation spectrum. We present here both the 
spectrum as a function of A, and as a function of k, see 
Fig.HUIandlTTl It differs drastically from that on Fig. [SJ 
Namely it developed an additional minimum at A sa 2Ao- 
We only plot numerical curves on these graphs as the an- 
alytic curves given by Eqs. (|7U|l . (j7T|l . and (J72J) are 
almost indistinguishable from the numerical ones. 

We can understand the appearance of a minimum at 
A = 2Ao by noticing that the particles which move at 
A = 2Ao scatter off the bulk of the Fermi sea exactly at 
resonance. Therefore, these particles try to form bound 
states with those in the Fermi sea and the resulting 
bosons then Bose condense. 

Now let us look the spectrum at Ao = 2.5. The density 
of states is again not much different from that on Fig. {7\ 
but the spectrum changes again, and is shown on Fig. 1121 
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FIG. 11: e(k) as a function of the momentum of the excita- 
tions k. Here K is taken from Eq. 1411 with ci = S2 = 0.01, 
and Ao = 3.3. 
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FIG. 12: e(A) as a function of A. Here K is taken from 
Eq. igTJ with ci = c 2 = 0.01, and A = 2.5. 



and Fig. ^] Remarkably, the spectrum develops areas 
where e < 0. By filling up the states where e is negative, 
it is possible to lower the energy of the system compared 
to that of the ground state. Therefore, the state dis- 
cussed in this section becomes unstable, and is no longer 
the ground state in this regime. It can be checked numer- 
ically that this instability develops roughly at Ao ~ 2.9, 
just as predicted analytically. This is why we refer to the 
state where only the states between — Xp and Xf arc filled 
as "simple ground state" . At sufficiently high density of 
the gas the ground state changes. 
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FIG. 13: e(k) as a function of the momentum of the excita- 
tions k. Here K is taken from Eq. I4H with ci = £2 = 0.01, 
and Ao = 2.5. 
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FIG. 14: The density of states p(A) for the case of large c, FIG. 16: l(k) as a function of the physical momentum k. 
Eq. 05). c 3 = 0.01, Ai = 0.65. Here K is taken from Eq. 05) with c 3 = 0.01, and Ai = 0.65. 
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FIG. 15: e(A) as a function of A. Here K is taken from _ 

Eq. 02) with c 3 = 0.01, and Ar = 0.65. FIG - 17: as a Action of A. Here K is taken from 

Eq. 05) with c 3 = 0.01, and Ai = 0.55. 



3. Numerical solution at large c 

If c is very large, then the kernel K can be taken from 
Eq. (|42J) . We are not able to find an analytical solution 
for this case, and present the numerical solution. It is 
qualitatively similar to the one at small c, although it 
has several new features. 

We take £3 = 0.01, and present numerical solution for 
Ai = 0.65, above the instability, and for Ai = 0.55, be- 
low the instability. Ai can generally take values smaller 
than appropriate values of Ao at small c because of the 
condition Eq. I|52|) as compared to Eq. (|51|) . The density 
of states at Ai = 0.65 is shown on Fig. Notice it is 
quite different from the analytically understood case of 
small c, depicted on Fig. 

As far as the energy spectrum, it is depicted on Fig. 1151 
It has minima similar to the ones on Fig. [5] In addition, 
not only particles, but also holes have nonmonotonous 
spectrum. Fig. 1161 shows the same spectrum, but as a 
function of a physical momentum k. 

Fig . 1 1 71 and 1181 show the excitation spectrum for Aq = 
0.55. Now the spectrum is unstable. 



VI. BEYOND THE SIMPLE GROUND STATE 

If the density of gas is increased beyond the critical 
density, at small c or at large c, then we can no longer 
use the Lieb-Liniger and excitation integral equations 
Eqs. 147|) and l|54|) . Instead, we can take a different path. 



The problem we confront is that above the critical den- 
sity we no longer know which states are filled and which 
states are empty in the ground states. We can circum- 
vent this problem by writing down the zero temperature 
version of the Yang- Yang equation ^l| , which looks like 

1 f°° A 2 
<A)-— J diiK{\-n)[e{ii)-\e{n)\] = --h. (78) 

This equation has to be solved at a given chemical po- 
tential h. Once its solution is known, we know that the 
particles states are those where e(A) > 0, while the hole 
states are where e(A) < 0. We can now determine the 




FIG. 18: l(k) as a function of the physical momentum k. 
Here K is taken from Eq. 1421 with S3 = 0.01, and Ai = 0.65. 
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density of states by solving 

PW ~ 4^ / d»K{\- [1 - sign e(A)] = i- . 

X (79) 
In case of low density resonant gas it is easy to solve 
Eq. H78|) as the integration there is effectively over the 
region where e(A) is negative. In those simple cases, we 
know that e(A) is negative only for — Xf < A < Xf and 
Eq. H78fl reduces to Eq. I|54|l . But above the critical den- 
sity, it is not at all clear how to solve these equations, ei- 
ther analytically or numerically. Numerically any naive 
scheme one could use to solve Eq. (|T8)) becomes unstable. 
For example, trying to find the solution by neglecting the 
integral in Eq. (|78|l first, and then taking it into account 
via iterations, immediately leads to an instability and a 
solution which quickly diverges with each successive it- 
eration. In fact, it is proven in Ref. ^l| that solving 
Eq. (|78J) by iterations works. However, their proof essen- 
tially uses the monotonicity of K(X) for A > 0, that is, it 
breaks down for resonant interactions. 

It is intuitively clear that above the critical density the 
single particle states whose momentum A falls into the 
areas of negative e(A), such as on Fig. El will be filled. 
But all of those states cannot be filled. Otherwise, the 
analog of the condition Eq. (|49|l (which is the very same 
condition, except the regions of integration are now over 
where e(A) is negative) would be violated. Unfortunately, 
it is not clear how to improve this intuition and make it 
more quantitative. 

Solving Eq. I|78|l and studying the resonant gas above 
the critical density will be left for future work. 



VII. CONCLUSIONS AND OUTLOOK 

We presented the technique for solving the problem of 
one dimensional resonantly interacting Bose gas exactly. 
With the Hamiltonian given by Eq. , supplemented by 
the condition Eq. J2J), we found the solution analytically 
for small c, or c <C \J eo/m a , at density below the critical 
density Eq. Q). In particular, we found the exact exci- 
tation spectrum, given by Eqs. |JH3J|, (7UJ), O, and G21- 
The excitation spectrum is not monotonous and has min- 
ima above the Fermi momentum. A typical excitations 
spectrum, as a function of asymptotic momentum A is 
shown on Fig. 1101 An excitation spectrum as a function 
of a physical momentum of a quasiparticlc is shown on 
Fig. ^] Above the critical density the solution we found 
becomes unstable and we were not able to find a stable 
solution in this regime. These features of the solution to 
the resonantly interacting gas remain valid at large c as 
well, as confirmed by numerical solutions of the Bethe 
Ansatz equations. 

The minima in the excitation spectrum should be eas- 
ily observable experimentally, if such a one dimensional 
gas is realized in a laboratory. Indeed, at any finite tem- 
perature, one should observe low energy excitations far 



above the Fermi momentum. These could be observed in 
a time of flight experiment. 

At c ~ y/ eo/m a , the excitation spectrum does not de- 
velop minima and the gas can be described by the Lut- 
tinger liquid theory. However, in one difference from the 
nonresonantly interacting Bose gas, its Fermi momentum 
Xf can never exceed a certain cutoff momentum in this 
regime. The increased density of gas is accommodated 
by the ever increasing density of states, as opposed to 
the Fermi momentum. 

In view of the possible future experiments, it is impor- 
tant to generalize the theory presented here to nonzero 
temperature. This would be straightforward to do by 
employing the appropriate Yang- Yang equation. 

The nature of the phase above the critical density for 
large or small c remains undetermined, and should be 
another subject of future work. 

A major limitation of the technique used in this paper 
is the neglect of three-body processes. Although justified 
in the relatively low density regimes considered here, a 
major technical improvement would involve coming up 
with a technique where this deficiency is absent. Per- 
haps an integrable version of Eq. Q can be formulated 
or some other technique can be proposed. Whether this 
is possible to do is currently not clear. In the absence of 
exact techniques, bosonization remains the only alterna- 
tive if one wants to study high density regimes where the 
density of molecules is supposed to be of the order of the 
density of atoms. 

Finally, an interesting generalization of the work pre- 
sented here would be to study a gas of fermions reso- 
nantly interacting and confined to one dimension. The 
main difficulty which will have to be overcome is that 
fermions carry spin, the 5-matrix becomes a real matrix 
in the spin space, and the Bethe Ansatz techniques be- 
come much more technically involved than in the case 
considered in this paper. 
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APPENDIX A: LOW ENERGY ONE 
DIMENSIONAL SCATTERING 

We would like to derive formulae Eq. (13} . Closely 
following Ref. ^2]> we observe that 

\S s , A \ 2 = \l + 2fs,A\ 2 = l- 
It follows from here that 
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and thus, Re f s A 



fs,. 



-1. We deduce 
1 

~ iF s , A (k)-l 



where Fs,a is a function of momentum k, real when k is 
real. Notice that 



cot 6s„ 



-F s , A {k). 



Next we would like to analyze the behavior of F at 
very small momenta. For that we need to study the 
Schrodinger equation at small energy. We observe that 
it is possible to split the range of x into several intervals. 
In the interval — tq < x < ro, we can neglect energy to 
find the Schrodinger equation 

1 d 2 ip 



+ U{x)ip = 0, |x| < r 



(Al) 



2m dx 2 

This equation must be solved with the boundary con- 
ditions ip(0) — 1, d ^f ^ — for symmetric amplitude, 



and with the conditions -0(0) = 0. 
antisymmetric amplitude. 



dx 



1 if we study 



In the interval rp 



< 



< 



1/k, we neglect the potential 



and the energy to arrive at 

1 d 2 ip 



= 0. 



2m dx 2 

Its solution is ij> = c\ + C2X. The constants c\ and C2 
can be determined by matching this solution with the 
solution to Eq. (|A1I) in the range of \x\ ~ tq. The cru- 
cial point is that they are energy independent. They of 
course depend on whether we are looking at symmetric 
or antisymmetric amplitude. 

Finally in the range \x\ > 1/k, the potential (but not 
the energy) can be neglected, and the Schrodinger equa- 
tion reduces to 

1 d 2 ip k 2 



2m dx 2 



2m 



The solution to this reads Acos(kx + 6s) in the sym- 
metric case or Asiri(kx + 6a) in the antisymmetric case, 



if 



> 



1/k and their symmetric (antisymmetric) reflec- 



tion if x < —1/k. Extending it to \x\ < 1/k, we find 
cot 6s = —kc\/c2- Thus Fs = kcl/c2 at small k, and we 
arrive at 

t.- 1 



ikF s (k 2 ) 



1 



where Fs(k) = Fs(k)/k, and Fs(0) is finite. Quite anal- 
ogously, taking the asymptotics Asui(kx + 6a), we find 
k cot 6a = c-i/c\, leading to 

fA 



~ k FA{k 2 



1 



where Fa(0) is finite. 

Finally, we observe that for purely imaginary k the 
scattering amplitude must be real. For that the functions 
Fs and Fa must be even functions of momentum k, which 
explains why they are the functions of k 2 . 



APPENDIX B: T-MATRIX CALCULATION OF 
THE SCATTERING AMPLITUDE 

We would like to calculate the T-matrix of the a- 
particles described by Eq. pD)l. The definition of T- 
matrix is given in many textbook, see, for example, 
Ref. 0| . It is basically the exact Green's function with 
free Green's function subtracted and external legs ampu- 
tated. 

We note that the free Green's function is given by 



G (E;k) = 



1 



and the Green's function of 6-particles is 

1 



D(E) 



E-eo + iO 

The T-matrix is given by the geometric series 
T = c + cGc + . 



1-Gc 

Here G represents trace of the Green's function 
roc dk 1 



G = 



and c is 



2tt e - I 1 + iO 

2m 



c = c + gDg. 



2E' 



Thus the T-matrix is 
c + gDg 



T = 



l-G(c + gDg) \c(E-e )+g 2 



E-e 



-G 



Upon substituting E = k 2 /2m and using Eq. I|21() we get 
Eq. 123- 



APPENDIX C: BETHE ANSATZ FOR THE 
MANY-BODY WAVE FUNCTION 

We would like to discuss why Eq. pHjl solves the ap- 
propriate Schrodinger equation. The proof is known from 
the studies of Bethe Ansatz. We are not going to present 
a mathematical proof, but merely list the reasons why 
this is the case. 

If x\ ^ x% 7^ ■ ■ • 7^ x n , then clearly this wave function 
solves the Schrodinger equation by construction, with the 
energy given by Eq. (|^T|) . 

Let us now take Xj close to Xk so that there are no 
other particles between these two. The terms in Eq. (|3Uf) 
can be split into pairs which differ by exchanging Xj and 
Xk- Let's take one such pair. Then if Xj < Xk, the wave 
function which includes the sum over this pair will be 
proportional to 



^ ^ e i\iXj+iX 2 x k _|_ e i\ix k +i\ 2 Xj g 



Ai — A2 
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where Ai and A2 are two momenta which Xj and Xk get 
multiplied by in the exponential in each of this pair of 
terms. If, however, Xj > Xk, then the wave function is 
proportional to 

Ai - A 2 \ e i\ 2 x 3 +i\ lXk 



Introducing relative coordinates 



R 



2 ' 



T — Xj Xk , 



Ai — A2 



A = Ai + A 2 



we find 



1> 



(e lAr + S{X)e' lXr ) 



for r < and 



AH 



— i\r 



s{xy 



for r > 0. This coincides with the definition of the 
scattering matrix Eq. (JSJ. So Eq. (|5U|l indeed solves 
the Schrodinger equation even when the particles pass 
through each other. 

Now let us show that if Xj = Xk = A, then ip, defined in 
Eq. (|29|) . is antisymmetric under the exchange of Xj and 
Xk- Indeed, suppose that Xj < Xk- Then once Xj and Xk 
are exchanged, we need to bring xt through all the par- 
ticles separating Xj and Xk until it occupies the position 
Xj normally occupies. At the same time Xj should be 
brought back to the position it normally occupies. Each 
exchange of Xk with the particles xi in between brings 
5 ( 2 ) ' wnere I labels these particles. Each exchange 
of Xj brings the factor 5 ( Ai 2~ A ) ■ Finally, an exchange of 
Xj and Xk brings the factor 5(0) . Thus ip, where Xj and 
Xk were exchanged, differs from ip before the exchange, 
by the factor of 



5(0) n s 



A -A, 



5 



A, -A 



-1. 



(CI) 



Here the product is over all the particles between Xj and 
Xk- And indeed, S(X)S(— A) = 1, as follows directly from 
the form 



5(A) = 



iXF(X 2 ) + 1 
iXF(X 2 )- 1' 



which follows from Eq. i|13|) . At the same time 5(0) = 
— 1, as follows from the same relation. 

This important property, from which the bosonic 
"Pauli principle" immediately follows, has one notable 
exception. Take the scattering amplitude Eq. 1(27(1 and 
set ei = 0. This corresponds to e — g 2 /c = 0, or in 
other word, to the situation where there is a bound state 
exactly at zero energy. Then, as can be easily checked, 
5(0) = 1. Thus in these circumstances there is no bosonic 
"Pauli principle" . Fortunately, this only happens when 
the condition of interest to us, Eq. @, is violated. 



APPENDIX D: GROUND STATE AND 
EXCITATIONS OF THE BOSE GAS WITH 
DELTA FUNCTION INTERACTIONS 

Consider the Bose gas Eq. (I28|) with repulsive two body 
potential U(x) = c 5(x) (c > 0), and with total density 
(number of particles per unit length) D. 

The solution to this problem can be found by solving 
the equations of Bethe Ansatz. Generally this involves 
solving integral equations which can only be solved nu- 
merically. The limits of very strong or very weak interac- 
tions can nevertheless be solved analytically. We present 
the solutions here. 

If the interaction is very strong (c ^S> D), then the 
kernel of the integral equation, Eq. H39|) vanishes. The 
Lieb-Liniger equation can be solved elementary, to give 

pW = ^- 
The total density of the gas is then 



D 



The Fermi velocity of the gas is given by (see Ref. [TU p 

D 



V F = -dD 



= X, 



and the excitation spectrum is given by 

e(A) = y-fc, 

where h is the chemical potential. 

Thus the Bose gas with strong repulsions is equivalent 
to a non-interacting Fermi gas, as is well known 

A somewhat more nontrivial regime is the one where 
the interaction strength is weak compared to density, 
c <C D. Then the Bethe Ansatz equations can be ap- 
proximately solved in the following way. 

First, consider the Lieb-Liniger equation Eq. 1)47(1 . with 
the kernel K given by Eq. It takes the form 

P(A) - f I" n 2CPi : ] 2 = f ■ 
2tt J_ Xf (A - vl) 2 + c 2 2tt 



First we rescale variables to find 



,MA) - -!- / dfi 



2cp(fi) 



1 



2tt 7_i (A - ft) 2 + c 2 2tt 

If c is very small, then we use an approximation 

1 25 c 

— -= w (5(A) + 

2vrA 2 + 5 2 ttA 2 



This gives 



dfi -p^Ei = A. 

-1 X-fl 2c 
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The solution to this well known integral equation states 
P(A) = ^^, (D2) 

or equivalently 



A 2 

A 2 



(D3) 



We emphasize that Eq. I|D3H is only correct for |A| < Xp, 
and it definitely breaks down for A close to the Fermi mo- 
mentum. For example, taken naively Eq. (|D3|I predicts 
that p(\f) — 0. In truth, p(Xf ) is not zero, but it is just 
much smaller than 1/c. 

Eq. (|D3|1 can be used to compute the total density 



r-A F 12 

D = I dfi p(fi) = 



4c 



From here we can compute the Fermi velocity 



D X F 
3X7 z 



(D4) 



(D5) 



Notice that this is exactly half of the Fermi velocity in 
the strong interaction case. We can now find the chemical 
potential h using 



dh 
dXp~ 



\ 2 



The excitation spectrum can be found by solving the in- 
tegral equation for the excitations 



(If I 



2ce(p) 



A 2 A 2 

2^ J-x F (A - m) 2 + c 2 = T ~ T' (D6) 



With the appropriate rescaling Eq. (|75|) . with the 
rescaled e(A) = q 2 e(X), and by taking the limit of small 
c, we bring this integral equation to the form 

1 X-fi 25 \ 6 4 J 
The solution to this equation reads 



e(A) 



A F 

Gc 



1 - 



, |A|<Ai 



(D7) 



For larger values of A we can use the integral equation 
Eq. (OBI to find 



e(A) 



A F A A 2 



T3--1, |A| >A F . 



(D8) 



Notice that in the vicinity of |A| = A^, both Eq. (|D7|) and 
Eq. (|D8|) break down. In particular, these expressions 
cannot be used to compute the Fermi velocity Eq. JD5I) 



simply by differentiating the energy e with respect to the 
momentum at A = A^?. 

It is also instructive to compute the physical momen- 
tum of the excitations fc, using Eq. (|55|l . For this purpose 
the phaseshift 9(X) can be replaced by 7rsignA, which 
gives 



k 



A 



r „ arcsin A . - , 
A + -Vl-A 2 + ^^, |A|<1, 

k = A+^, |A| > 1. (D9) 



Here k = 1c/Xf- 

Finally, the Luttinger parameter of the Bose gas is 
given by 



9 



ttD 

VF 



(D10) 



The second equality is valid only in the limit j » 1. 
The Luttinger parameter at small c can also be derived 
with the bosonization techniques |l6j| . and the answer 
coincides with Eq. i|D10|) . 



APPENDIX E: MOLECULAR DENSITY AND 
THREE-BODY PROCESSES 



We would like to give an estimate of what percent- 
age of the atomic wave function is concentrated inside 
the molecules. We will do it by directly calculating the 
wave function of the single-particle Fano Anderson model 
Eq. (|20|l . Indeed, its Hamiltonian is quadratic, so it can 
be diagonalized by the appropriate rotation of the oper- 
ators a, b. Its eigenstates at the energy e can be con- 
structed by writing 

V(e) = c 6+ |0) + £> fc <4 |0) . 

k 

where Co and ctk satisfy, curtesy of Eq. (|20l) , the following 
equations (L is the size of the system) 



eo c 
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E 



ctk = e c , 



'VI 

It follows that 



9 c \ - 



c e — ei 
a a = — — c , 



g\fZ e - e q 
while the energy level e satisfies 
c \ -> e — c\ 



= e - e . 



(El) 



Here t\ is given, as before, by Eq. I|25[) . and e q = q 2 / (2m). 
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The normalization condition for the wave function 
states 



9 




Thus the condition that the particle spends most of its 
time outside the molecular state, created by 6 T , can be 
formulated as follows 



— Y 



> 1. 



(E2) 



We would like to apply this condition to make sure that 
atoms spend only a small fraction of their time bound 
into molecules in the regimes of interest in this paper. 

Consider first the case of small c. At densities up to 
the critical density Eq. the particles from the Fermi 
sea, having momenta from — Xf to A_f, are not in res- 
onance with each other, so their quasistationary bound 
states can be completely neglected. Suppose we add one 
more particle. Its spectrum is shown on Fig. 1101 If its 
momentum is of the order of 2Ao, the particle is in reso- 
nance with particles in the Fermi sea and spends at least 
some of its time bound into a molecule. This corresponds 
to the local minimum in e(A) as shown on Fig. 1101 In or- 
der to be able to trust the physics of Fig. we need to 
verify that the condition Eq. i|E2(l holds for the major- 
ity of the two body processes between the particle with 
momentum close to 2Ao and the particles in the Fermi 
sea. 

We observe that at c small, essentially only one term 



with optimal q, where e — e q 



contributes to the sum 



in Eq. QElj l. We can use that to estimate e — e q and 



substitute into Eq. i|E2|) . again retaining only one term 
in the sum. We find 



L ( 



e y » i. 



This condition is clearly violated for the particles exactly 
at resonance, so that e = eo- Yet we only need to check 
that this condition holds for the typical pair of particles, 
not just those which are exactly at resonance. The typical 
momentum difference for a particle at A = 2Ao and a 
particle in the Fermi see is given by 2Ao + Xf, and a 
typical energy is 



e - £o = 



(A 



2m 



\ 2 
A o 

2m 



Here we used the estimate Eq. (jHJ for Xf- Hence the 
criterion becomes 

L 2 s i 

— C3£,?m3 > 1. 
9 

Observing that g 2 ~ c (due to ei, defined in Eq. 
being finite) we find that the above condition definitely 
holds at small c. Thus we proved that at small c, in the 
regime of interest to us where the density is of the order 
of Eq. the particle spends the overwhelming majority 
of its time outside the molecular state. 

The case of large c is somewhat harder to analyze as we 
do not have an analytic solution of this problem even with 
two-body interactions neglected. Nevertheless we expect 
the condition Eq. (|E2I) to hold as c is now large, and, 
unless the energy e is exactly at resonance, the particle 
still spends most of its time outside the molecular state. 
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